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Abstract. Bose-Einstein correlations of like-sign charged hadrons produced in deep-inelastic electron and 
positron scattering are studied in the HERMES experiment using nuclear targets of 1 H, 2 H, 3 He, 4 He, N, 
Ne, Kr, and Xe. A Gaussian approach is used to parametrize a two-particle correlation function determined 
from events with at least two charged hadrons of the same sign charge. This correlation function is compared 
to two different empirical distributions that do not include the Bose-Einstein correlations. One distribution 
is derived from unlike-sign hadron pairs, and the second is derived from mixing like-sign pairs from different 
events. The extraction procedure used simulations incorporating the experimental setup in order to correct 
the results for spectrometer acceptance effects, and was tested using the distribution of unlike-sign hadron 
pairs. Clear signals of Bose-Einstein correlations for all target nuclei without a significant variation with 
the nuclear target mass are found. Also, no evidence for a dependence on the invariant mass W of the 
photon-nucleon system is found when the results are compared to those of previous experiments. 
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Introduction 

Hadron production in deep-inclastic scattering (DIS) of 
leptons off nuclei is a powerful tool to study the quark 
hadronization process. The distance scale over which a 
struck quark that received a sufficiently large energy- 
momentum transfer from an incident lepton develops into 
a colorless hadronic particle extends well beyond the size 
of a single nucleon. Therefore, the distribution of hadrons 
in the final state may be modified by interactions of the 
developing hadronic state with the nuclear medium out¬ 
side the struck nucleon. In general, this intermediate state 
is some mixture of quarks and gluonic fields that have 
not reached their asymptotic (confined) states, and so 
any modification should depend on the evolution of that 
state. Similarly the fully formed hadron may still pass 
through the nuclear medium and be subject to rescatter¬ 
ing processes (see, e.g., Ref. [1]). 

One means of studying the final hadronic state is 
the use of Bose-Einstein correlations (BEC) in the dis¬ 
tribution of bosons, in particular pions. These correla¬ 
tions arise from interference between different parts of 
the symmetrized wave function of identical bosons from 
incoherent sources. This well-known technique of inten¬ 
sity interferometry was first developed by Hanbury Brown 
and Twiss to measure stellar radii [2]. Its first use in 
particle physics, half a century ago, was to study the 
pp annihilation process [3,4] with incident anti-protons 
of 1 GeV momentum. Since then many measurements 
of BEC have been performed in hadron-hadron scatter¬ 
ing experiments. In addition, several studies of BEC in 
the e + e - annihilation process have been performed (see, 
e.g., Ref. [5]), especially by the LEP experiments. Mea¬ 
surements of BEC from deep-inelastic lepton scattering 
experiments are less abundant. The results from experi¬ 
ments using charged leptons as incident particles can be 
found in Refs. [6,7,8,9,10], while the results from neu¬ 
trino experiments are found in Refs. [11,12,13]. Several 
reviews [5,14,15,16] summarize the present theoretical 
and experimental knowledge of BEC. The theory of BEC 
in particle physics was originally developed in the papers 
of Kopylov and Podgoretskii [17,18,19] and Cocconi [20]. 
It should be noted that most of the theoretical work has 
focused on the understanding of BEC in heavy-ion col¬ 
lisions, in which a “fireball” source distribution, created 
by the collision roughly at rest involving many parton el¬ 
ementary interactions, decays into hadrons. Only a few 



Fig. 1. Schematic illustration of the Bose-Einstein effect. 

references consider the quite different case of fragmenta¬ 
tion in DIS and e + e - processes, in which quite different 
hadron-momentum and spatial-source distribution might 
be assumed (see, e.g., Ref. [21,22]). Estimates of BEC in 
e + e - annihilation from string-fragmentation models [22] 
indicate that correlation parameters are mostly depen¬ 
dent on string-breaking parameters, because the strongest 
correlations are from pions resulting from adjacent breaks 
along a string. 

To better understand the underlying physics of BEC 
one may consider a simple example of the emission and 
detection of two identical bosons, e.g., two pions, from 
points r a and rp, which are observed with momenta k a 
and kb at detectors a and b (Fig. 1). The two pions are 
indistinguishable and the total wave function of the two- 
pion system must be symmetric under the exchange of 
them: 

KrMli + %c&apj , ( 1 ) 

where \^ aa is the wave function of a pion produced at 
point r a and observed at detector a while is the 
wave function of a pion produced at point rp and ob¬ 
served at detector b. Assuming plane waves, i.e., tPaa ~ 
exp (ik a r a ), one may obtain Itf^l 2 = 1+cos (Sk-Sr) with 
Sk = k a — kb and Sr = r a — rp. Thus the correlation 
function resulting from the interference of the two terms 
in Eq. (1) will take the following form: 

R(k ai kb) oc 1 + cos(6k ■ dr). (2) 

This expression shows that the BEC effect measures the 
projection of the spacial distance (Sr) between two par¬ 
ticle sources on the direction of the momentum difference 
(Sk) between the observed pions. One can generalize two- 
point sources to a continuous space-time distribution of 
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sources. Experimentally this is achieved by finding the 
boson correlation function. In BEC this correlation ob¬ 
servable is defined in terms of the two-particle correlation 
function 

R(PuP 2) = D(p 1 ,p 2 )/[D(p 1 ) ■ D(p 2 )} , (3) 

where pi and p 2 are the particle four-momenta, D(pi,p 2 ) 
is the two-particle probability density and D(pi), D(p 2 ) 
are one-particle probability densities. 

Typical analyses use models for the correlation func¬ 
tion R with a limited set of parameters. In the case of fire¬ 
ball decay (heavy-ion collision), the correlation function 
might be described by a (spherical) Gaussian distribution 
in space with an independent parametrized exponential 
decay in time. However, such an approach is not Lorentz 
invariant. Generally the Goldhaber parametrization [4] is 
a more convenient approach, in which the distribution 
is a function of the Lorentz-invariant quantity T 2 of the 
hadron pair, where T 2 = — (pi — p 2 ) 2 = S — 4m 2 with S 
the squared invariant mass of the pair and m n the pion 
mass. In this parametrization the correlation function R 
takes the form 

R(T) = l + A-e" r2 ^. (4) 

This parametrization corresponds to a Gaussian shape of 
the particle source distribution of size re in the center of 
mass of the pair. The additional parameter A, the chaotic- 
ity (or incoherence), has been introduced to account for 
a possible incoherent contribution of the pion emitters. 
Completely coherent pion sources lead to the absence of 
correlations among the bosons (A = 0) [23,24], while in 
the simplest theoretical treatment, completely incoherent 
sources lead to A = 1. Subsequent theoretical work has 
shown that a number of different effects can modify the 
value of A in both directions, greater and smaller than 
unity (see, e.g. Ref. [14]). Examples of such effects are 
the presence of decay products of long-lived resonances 
in addition to directly produced bosons in the final state, 
final-state interactions (FSI), and the deviation of the 
exact shape of the source distribution from the assumed 
Gaussian form. Other experimental effects can also influ¬ 
ence the experimentally measured value of A, often due 
to the purity of the boson sample as determined by the 
quality of the particle identification within the experi¬ 
ment. For these reasons, it is difficult to compare the re¬ 
sults of A between different experiments, even at identical 
kinematic conditions. 

Note that the Goldhaber parametrization does not 
correspond to any well developed theoretical approach 
but rather serves as a convenient tool for the comparison 
of experimental results. The parametrization describes 
the shape of the distribution in T satisfactorily. How¬ 
ever, depending on the physics of the reaction, a rigor¬ 
ous treatment of the coherent and incoherent sources re¬ 
quires modification of the functional form of the correla¬ 
tion function (see, e.g., Ref. [25]). 

Given these remarks, one must note that there are 
additional points to consider in a full treatment of BEC. 
Initial correlations among the bosons are affected by FSI 


between the produced particles as well as with the pro¬ 
duction environment, both via the Coulomb interaction 
and hadronic interactions. The long-range Coulomb FSI 
is quite often accounted for by introducing a multiplica¬ 
tive Gamow factor [23]. It increases (decreases) the two- 
particle density for opposite-(like-)sign particles. The cor¬ 
rection factor is essential only at very small values of T 
and very quickly approaches unity with increasing T. For 
pion pairs at T = 0.05 GeV the correction factor differs 
from unity by only about 1.5%. There are other more 
elaborate calculations that predict an even smaller mag¬ 
nitude for the required correction and also examine its 
model dependence. Short-range strong FSI between the 
two identical pions may also influence their correlation 
(see, e.g., Ref. [5,14,15,16,26]); without clear theoretical 
estimates for the kinematics of the present experiment it 
was chosen not to attempt any correction for both long- 
range and short-range FSI. 

Effects that can alter the pair correlation within the 
nuclear environment are the focus of this study of BEC 
for nuclear targets ranging from hydrogen to xenon. Mea¬ 
surements with the same experimental apparatus and 
kinematics help to minimize possible systematic bias in 
the observed target-mass dependence. In DIS a difference 
in the size of the particle emission region could exist for 
pions produced off a free nucleon as compared to that off 
bound nucleons. In addition, interactions of the struck 
bare quark, during fragmentation, or of the fully devel¬ 
oped hadron with the nuclear environment could alter 
both the apparent size of the emission region and the 
amount of incoherence. For example, a common simple 
assumption is that the correlations of a pair of identical 
pions are determined by the relative positions of their last 
scattering points, which thereby play the role of indepen¬ 
dent particle sources. If the pions scatter from the nuclear 
matter, one would expect an increase in the size of the 
emission source as a function of target radius. Another 
example would be the increasing probability of gluon ra¬ 
diation within the nucleus leading to a change in the pion- 
pair correlations relative to hadronization in free space. 

The influence of nuclear re-scattering processes on 
BECs in heavy-ion collisions was studied in Ref. [27] and 
effects on the source size of 15-20% were found. No es¬ 
timates exist for the case of lepton-nucleus hadron pro¬ 
duction. Earlier experimental studies of BECs from DIS 
by nuclei are quite limited. The BBCN Collaboration [12] 
found the BEC parameters tq and A to be independent 
of the atomic mass. These measurements, however, are 
limited to three light nuclei, 1 H, 2 H, and Ne. 

There is significant evidence that the nuclear medium 
affects the hadronization processes in lepton scattering 
at HERMES energies. Results on this subject have been 
presented in a number of papers [1,28,29,30,31,32,33, 
34]. In particular, the transparency of the 2 H, 3 He, and 
14 N nuclei to exclusive incoherent p° electro-production 
has been measured [28] and significant dependences on 
the coherence length were found. A series of papers [1, 
29,30,31] is devoted to the investigation of the hadron 
multiplicity variation in different kinematic regions and 
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its dependence on the target atomic mass A up to xenon. 
The most prominent features of the data are an increased 
hadron attenuation with increasing value of the mass 
number A of the nucleus and the attenuation becoming 
smaller (larger) with increasing values of v ( z ), where v is 
the energy of the virtual photon in the laboratory system, 
and z = Eh/v is the fractional hadron energy. 

The influence of the nuclear medium on the ratio of 
double-hadron to single-hadron yields in DIS was also in¬ 
vestigated [32]. Nuclear effects are clearly observed but 
with substantially smaller magnitude as well as reduced 
A dependence compared to the single-hadron multiplicity 
ratios. The first detailed study of the dependence on the 
target nuclear mass of the average squared transverse mo¬ 
mentum (p 2 ) of hadrons produced in deep-inelastic lepton 
scattering is described in [33]. It is found that the aver¬ 
age squared transverse momentum is increasing with the 
atomic mass number. 

In short, several studies with different nuclear targets 
performed at HERMES show significant modifications of 
the hadron observables within the nuclear medium as 
compared to the results on a proton/deuteron target. 

Experiment 

The present measurement of the BEC was performed 
with the HERMES spectrometer [35] using the 27.6 GeV 
polarized lepton (electron/positron) beam stored in the 
HERA ring at DESY. The spectrometer consisted of two 
identical halves above and below the lepton beam line. 
The scattered lepton and the produced hadrons were de¬ 
tected within an angular acceptance of ±170 mrad hori¬ 
zontally, and ±(40-140) mrad vertically. 

All the targets were internal to the lepton storage ring 
and consisted of polarized or unpolarized 1 H, 2 H, and 
3 He, or unpolarized 4 He, N, Ne, Kr, or Xe gas injected 
into a thin-walled open-ended tubular storage cell. Target 
areal densities up to 1.4 x 10 16 nucleons/cm 2 were ob¬ 
tained for unpolarized gas corresponding to luminosities 
up to 3 x 10 33 cm -2 s _1 . The luminosity was measured 
using elastic scattering of the beam leptons off the elec¬ 
trons in the target gas, Bhabha scattering for a positron 
beam and Mpller scattering for an electron beam [36]. 

The trigger was formed by a coincidence between the 
signals from three scintillator hodoscope planes, and a 
lead-glass calorimeter where a minimum energy deposit 
of 3.5 GeV (1.4 GeV) for unpolarized (polarized) target 
was required. The scattered leptons were identified using 
a transition-radiation detector, a scintillator pre-shower 
counter, an electromagnetic calorimeter, and a threshold 
gas Cherenkov counter. In 1998 the threshold Cherenkov 
counter was replaced by a ring-imaging Cherenkov detec¬ 
tor (RICH). 

Analysis 

Scattered leptons are selected by imposing constraints 
on the squared four-momentum of the virtual photon, 


Table 1. Number of DIS events with more than one detected 
hadron, N ev , the number of like-sign hadron pairs, and 

of unlike-sign hadron pairs, N unl ' lke j that meet the kinematic 
requirements for each target. 


Nucleus 

N ev 

jylike 

jyunlike 

4 H 

1145046 

478946 

958185 

2 H 

1297356 

680143 

1178797 

3 He 

34391 

15295 

29165 

4 He 

79776 

30539 

59244 

N 

92968 

41112 

78402 

Ne 

175594 

75898 

146145 

Kr 

211456 

91391 

172946 

Xe 

106274 

46130 

87125 


Q 2 > 1 GeV 2 , and on the invariant mass of the photon- 
nucleon system, W 2 > 10 GeV 2 . The constraint on W 2 is 
applied in order to ensure the predominance of multiple 
particle production in the DIS events. 

The results of this study are based on data collected 
by the HERMES Collaboration between the years 1996 
and 2006. The yields from polarized targets are summed 
over both target spin orientations. The yields from all 
targets are summed over both (longitudinal) beam polar¬ 
ization states. 

Events with only one identified lepton of the same 
charge as the beam lepton and momentum larger than 
3.5 GeV are accepted. The presence of at least two charged 
hadrons with momenta 2.0 GeV < Ph < 15 GeV is re¬ 
quired for further analysis. The total numbers of such 
multi-hadron DIS events, N ev , the numbers of like-sign, 
N hke , and unlike-sign, N unhke ; hadron pairs available for 
the analysis are given in Table 1 for each target. 

In this analysis all charged hadrons are considered to 
be pions. Simulations using the PYTHIA 6.2 event gener¬ 
ator [37] , tuned to provide an accurate description [38] of 
semi-inclusive deep-inelastic hadron lepto-production in 
the HERMES kinematic region, show that the observed 
charged hadrons are distributed in relative proportion 
tt/K/p(p) = 78%/12%/10%. In the same simulations one 
finds that 55% of like-sign hadron pairs and 66% of unlike- 
sign hadron pairs are truly pion pairs. This results in a 
“dilution” of the parameter A under the assumption that 
the non-pion pairs do not contribute to the BEC. The 
value of A in this analysis is expected to be smaller than 
about 0.5. Kaon like-sign pairs contribute only 2% while 
their unlike-sign pairs contribute about 4%. 

Experimentally, it is difficult to measure the inclu¬ 
sive single particle spectrum required to determine the 
probability density D(p) for all possible momenta p in 
the formal definition of the correlation function R(jp\,p 2 ). 
A common practice is to substitute the two one-particle 
probability distributions D(j>i)-D(j> 2 ) with a two-particle 
probability density reference distribution D r (pi ,P 2 )- This 
reference distribution is constructed from experimental 
two-particle distributions that do not have any BECs. 
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The experimental correlation function is then defined as 

R(Pi,P2) = D{p 1 ,p 2 )/D r {p 1 ,p 2 ). (5) 

The Goldhaber parametrization usually used includes an 
additional normalization parameter 7 and a polynomial 
function V(T) to describe the long-range correlations at 
large T and has the form 

R(T) = 1 -[l + \-e- T2r e)-V(T). ( 6 ) 

The long-range correlations at large T may appear due to 
charge and energy conservation, phase-space constraints 
and imperfections in the reference sample. The form of 
the polynomial V{T) used to model these long-range cor¬ 
relation effects is taken to be (1 + 8 ■ T 2 ). A linear depen¬ 
dence (1 + S' ■ T) is used to estimate the influence of the 
chosen form on the final results. The parameter 6 and S' 
are free parameters like re, A, and 7 . 

The magnitude of the two-particle BEC is measured 
by comparison of the experimental distribution with a 
reference sample distribution [see Eq. (5)]. The method 
of constructing the reference sample is the main source 
of systematic uncertainty, especially since the multiplicity 
of hadrons in the HERMES experiment is relatively low. 
One of the main problems in measuring BEC is the eval¬ 
uation of biases caused by an imperfect reference sample, 
thus it is desirable to use at least two different approaches 
to construct a reference sample to cross-check the corre¬ 
lation results. 

Two of the most widely used methods to construct a 
reference sample are employed here: 

— Method of event mixing {MEM), 

— Method of unlike-sign pairs {MUS). 

In MEM hadron-pair distributions of the same charge 
are created by using hadrons from different events, while 
in MUS hadrons with different charge from the same 
event are used. Other methods can be found in the litera¬ 
ture, each with its own shortcomings. The main technical 
difficulty of the two methods chosen here is the violation 
of momentum and energy conservation in the kinematic 
event topology when selecting two hadrons from differ¬ 
ent events in the case of MEM , and the contribution 
of events from resonances that are not present in the 
like-sign distribution in the case of MUS. The system¬ 
atic effects associated with these two reference samples 
are studied using the PYTHIA-based Monte Carlo simu¬ 
lation of the HERMES experiment discussed above, the 
inclusion of quantum interference of BEC not being en¬ 
abled in the PYTHIA event generation. 

The construction of the reference sample using the 
MUS is done by forming the distribution in T of all 
unlike-sign hadron pairs, requiring the same constraints 
as for the like-sign pairs. For the MEM , to construct 
a sample of uncorrelated hadron pairs, a combination of 
charged hadrons from two different DIS events is used. 
The first hadron of a pair is taken from one event while 
the second hadron is taken from another event. Care must 
be taken to conserve collinearity of the virtual-photon 



Fig. 2. Top panel: normalized experimental (stars) and sim¬ 
ulated (line) distributions for unlike-sign hadron pairs as a 
function of the variable T. Bottom panel: the ratio of these 
experimental (exp) and simulated (MC) distributions. 


vectors <71 and q 2 from these two different events. For 
each event, the momentum vector of the total hadronic 
system must lie in the direction of the virtual photon 
q = p e — p' e , where p e and p' e are the momenta of the 
incident and scattered leptons. To conserve collinearity 
the momenta of all hadrons in the second event must be 
rotated in such a way that the total hadronic momentum 
is aligned along the direction of q± of the first event. 

The quality of the HERMES simulation with respect 
to the description of the measured unlike-sign hadron pair 
distribution is demonstrated in Fig. 2 for a sample of 
hydrogen target data. The top panel shows the experi¬ 
mental T distribution of unlike-sign hadron pairs (exp) 
in comparison with the simulated data for h + h~ pairs 
(MC). The bottom panel shows their ratio. The figure 
demonstrates good agreement between the experimental 
and simulated distributions. Based on this agreement, the 
simulation results are used to further reduce the system¬ 
atic biases of the reference sample and experimental dis¬ 
tributions through the use of a double-ratio definition for 
the correlation function R{T). For the two methods in 
this analysis, 

rMEM _ {Hke/ r ni Xe d) exp / {like/mixed) MC , 

rMUS _ {Uke/unlike) exp / {like/unlike) MC . ^ ^ 

Dividing the experimental ratios by the simulated ratios 
is expected to reduce biases resulting from the violation of 
kinematic constraints in the MEM and from resonance 
contamination in the MUS, since these biases also exist 
in the simulated event distributions. 
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Fig. 3. Consistency check of the two chosen reference samples. 
The quantity R TST is defined in the text. The curve is a linear 
fit to the data for T between 0.05 GeV and 1.3 GeV. 

As a test of the validity of R(T) using the simulation 
results with a double ratio, the MEM was used with 
unlike-sign hadron pairs from the hydrogen data sample 
to construct the double ratio 

E^ST = (unlike / mixed) exp / (unlike/mixed) MG , (8) 

shown as a function of T in Fig. 3. This test ratio is 
expected to have no BECs, and ideally would have a 
value of unity over the entire T range. At very low T 
(< 0.05 GeV), at T « 0.4 GeV, and at T > 0.9 GeV 
this double ratio deviates from unity significantly. As 
shown by a linear fit to R TST there is a slight linear 
dependence over most of the range of T, indicating some 
small residual bias. The deviation near 0.4 GeV in the 
simulation is likely due to insufficient description of K$ 
production, which contributes to the N unllke distribu¬ 
tions (see Fig. 2). The deviations at very low and at 
large T likely arise from some combination of effects in 
both the simulation of the MUS and the MEM con¬ 
struction of the reference sample. The very low T region, 
T < 0.05 GeV, of the double ratio distributions is ex¬ 
cluded from further analysis due to lack of statistics. A 
fit to the correlation function R TST (shown in Fig. 3) 
with the Goldhaber parametrization [Eq. (6)] over the 
range 0.05 GeV< T < 1.30 GeV gives A=0.000±0.003 
and tg =0.0±1.4 fm, suggesting that the fluctuations at 
large T and the slight non-zero linear dependence on T do 
not cause a significant bias of the extracted parameters 
A and tq- 

Results 

The double-ratio correlation functions obtained from hy¬ 
drogen data are shown in Fig. 4 for both types of the 
reference sample. The curves in the figure are results of 
fits using the Goldhaber parametrization [Eq. (6)]. The 



Fig. 4. Double ratio correlation function for like-sign hadron 
pairs obtained with MEM and MUS based on hydrogen tar¬ 
get data. 


fits are performed over the range of 0.05 GeV < T < 
1.30 GeV. The values for the two parameters obtained 
from the fits are given in Table 2. 

The systematic uncertainties are estimated by varia¬ 
tions of the fit range in T, the bin width, and the polyno¬ 
mial form for the long-range correlations term, i.e., using 
a linear dependence (1 + S'T). The results of the two 
different methods are consistent (see Table 2). Values of 
the fit parameter 5 from the quadratic form of V(T) are 
—0.08 ±0.01 and —0.05 ±0.01 respectively for the MEM 
and MUS. 

The kinematic dependence of the BEC parameters 
on the invariant mass W of the photon-nucleon system 
has been studied for the hydrogen target data sample. In 
Fig. 5 the resulting parameters ro and A are presented 
for like-sign hadron pairs as a function of W obtained 
with the MEM and MUS methods. Within the present 
systematic and statistical uncertainties there is no clear 
dependence of the parameters on the invariant mass W 
in this range. Previous measurements from the HERA HI 
experiment [8] over a broad range at high W (65 GeV < 


Table 2. Results for the Goldhaber parametrization fitted 
to the HERMES hydrogen data, both for the mixed-event 
method (MEM) and the method of unlike-sign pairs ( MUS ). 


Method 

Goldhaber parameters 

MEM 

tg = 0.64 ± 0.03(stat)!g;o4(sys) fm 

A = 0.28 ± O.Ol(stat)Igjgs (sys) 

MUS 

tg = 0.72 ± 0.04(stat)lQ;og(sys) fm 

A = 0.28 ± 0.02(stat)lg Q 4 (sys) 
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W [GeV] 


Fig. 5. Parameter r G (top panel) and A (bottom panel) as a 
function of W, obtained with MEM and MUS methods on 
hydrogen. The inner and outer error bars indicate the statis¬ 
tical and total uncertainties. For the latter the statistical and 
systematic uncertainties are added in quadrature. 


W < 240 GeV) found only slight evidence of an increase 
in r G . 

As mentioned above, BEC has been studied in a num¬ 
ber of lepton-hadron and e + e _ experiments. The Gold- 
haber parametrization is used in most of these analyses. 
The parameter r G as a function of the average value of 
W in lepton-nucleon scattering experiments is shown in 
Fig. 6. The parameter A for a given experiment may de¬ 
pend on the hadron fractions and on the experimental 
details, hence the results of A obtained here are not com¬ 
pared to those in other measurements. In the majority of 
these experiments the extracted values of r G depend upon 
the method of the construction of the reference sample. 
Even for a single experiment, e.g., EMC, the parame¬ 
ter r G obtained with the MUS is twice as large as that 
obtained with MEM. From Fig. 6 no clear dependence 
of the parameter r G on W can be deduced, from neither 
methods (MEM and MUS). The following conclusions are 
drawn from a comparison of these results from the differ¬ 
ent experiments: 

1 . Most values of the parameter r G are in the range of 
0.4 fm to 1.0 fm. 

2. The results strongly depend on the choice of the refer¬ 
ence sample. Analyses of the same data set with differ¬ 
ent reference samples often give incompatible results 
for r G (and A). 

3. The MUS typically gives higher values for the param¬ 
eter r G than the MEM. 

The HERMES results on hydrogen are in general agree¬ 
ment with those of previous lepton-nucleon scattering ex¬ 
periments over a broad range in W, and agree well with 
the BBCN neutrino experiment, which is at a slightly 


2 
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Fig. 6. Goldhaber radius r G , as a function of W, obtained 
in lepton nucleon scattering experiments [6,7,8,9,11,12,13]. 
Different markers are used to indicate the different methods 
for the construction of the reference sample and the kinds of 
uncertainties included. 


higher mean W than HERMES. Similar results are seen 
in e + e _ collisions at LEP (see Ref. [5]). 

A possible nuclear dependence in BEC was examined 
using an extensive HERMES data set (cf. Table 1). The 
correlation function for like-sign hadron pairs produced 
in scattering off the nuclear targets 2 H, 3 He, 4 He, N, Ne, 
Kr, and Xe was determined using the same approximate 
parametrization as given in Eq. 6. Systematic uncertain¬ 
ties are estimated separately for each target and each 
reference sample (MEM and MUS). The parameters r G 
and A are presented in Fig. 7 as a function of the tar¬ 
get atomic mass A. No dependence of these parameters 
on target atomic mass is observed within the estimated 
uncertainties. Fit results with a constant over the whole 
range of the atomic mass for the four sets of data points 
are presented in Table 3. Here, the total uncertainty of 
each particular point is taken as the quadratic sum of sta¬ 
tistical and systematic uncertainties. The parameters ex¬ 
tracted with the two reference samples are in good agree¬ 
ment. 


Table 3. Fit of a constant to the Goldhaber parameters as 
a function of the target atomic mass A. Results are given for 
both the mixed-event method (MEM) and the method of 
unlike-sign pairs (MUS). 


Method 

Value 

X 2 / NDF 

MEM 

r G = 0.634 ± 0.017 fm 

1.5 


A = 0.289 ± 0.006 

2.1 

MUS 

r G = 0.636 ± 0.021 fm 

1.2 


A = 0.289 ±0.011 

1.4 














Fig. 7. The parameters tg (top panel) and A (bottom panel) 
are shown as a function of the target atomic mass A. The inner 
part of the error bars indicate the statistical uncertainty and 
the total error bars have systematic uncertainties added in 
quadrature. The horizontal lines correspond to the average 
value of the parameters. 


To date there are no theoretical estimates for the mag¬ 
nitude of nuclear effects on BEC in DIS. In the absence 
of some hitherto unknown effect of multi-particle corre¬ 
lations, hadrons produced are expected to interact with 
the nuclear medium. Within the sensitivity of this exper¬ 
iment no clear dependence of the parameters A and re on 
the target atomic mass is observed, consistent with earlier 
results by the BBCN Collaboration [12]. This is similar to 
the rather weak dependence of the double-hadron yields 
on the target atomic mass observed at HERMES [32], in 
contrast to much stronger effects observed in the distri¬ 
butions of single-hadron yields [1,29,30,31]. 

In conclusion, a study of the Bose-Einstein correla¬ 
tions between two like-sign hadrons produced in semi- 
inclusive deep-inelastic electron/positron scattering off 
nuclear targets ranging from hydrogen to xenon has been 
carried out. Two different methods of constructing the 
reference sample are used in this study, and Bose-Einstein 
correlations are clearly observed in all the data samples. 
The results obtained using the two reference sample meth¬ 
ods are in good agreement, suggesting that most of the 
systematic uncertainties connected with the construction 
of the reference samples are taken into account by the use 
of double ratios corrected via an accurate experimental 
simulation. Within the total experimental uncertainties, 
no dependence of the parameters re and A on the target 
atomic mass is observed. 
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